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Abstract. We present a study of the baryonic fraction in galaxy clusters aimed at constraining the cosmological parameters 
Clm, Oa and the ratio between the pressure and density of the "dark" energy, w. We use results on the gravitating mass profiles 
of a sample of nearby galaxy clusters observed with the BeppoSAX X-ray satellite (Ettori, De Grandi & Molendi, 2002) to 
set constraints on the dynamical estimate of Q m . We then analyze Chandra observations of a sample of eight distant clusters 
with redshift in the range 0.72 and 1.27 and evaluate the geometrical limits on the cosmological parameters f2 m , Qa and w 
by requiring that the gas fraction remains constant with respect to the look-back time. By combining these two independent 
probability distributions and using a priori distributions on both fib and Ho peaked around primordial nucleosynthesis and 
HST-Key Project results respectively, we obtain that, at 95.4 per cent level of confidence, (i) w < —0.49, (ii) fi m = 0.34^q' ( 1 )5, 
S7a = 1.30jl5'o9 for w = — 1 (corresponding to the case for a cosmological constant), and (iii) f2 m = 1 — Q,a = 0.331q 'q 7 5 
for a flat Universe. These results are in excellent agreement with the cosmic concordance scenario which combines constraints 
from the power spectrum of the Cosmic Microwave Background, the galaxy and cluster distribution, the evolution of the X-ray 
properties of galaxy clusters and the magnitude-redshift relation for distant type la supernovae. By combining our results with 
the latter method we further constrain Qa = 0.94lto §1 an d w < —0.89 at the 2a level. 

Key words, galaxies: cluster: general - galaxies: fundamental parameters - intergalactic medium - X-ray: galaxies - cosmol- 
ogy: observations - dark matter. 



1. INTRODUCTION 

Several tests have been suggested to constrain the geometry 
and the relative amounts of the matter and energy constituents 
of the Universe (see recent review in Peebles & Ratra 2002 
and references therein). A method that is robust and comple- 
mentary to the others is obtained using the gas mass fraction, 
/gas = M gas /M to t, as inferred from X-ray observations of 
clusters of galaxies. In this work, we consider two independent 
methods for our purpose: (i) we compare the relative amount of 
baryons with respect to the total mass observed in galaxy clus- 
ters to the cosmic baryon fraction to provide a direct constraint 
on Q m (this method was originally adopted to show the crisis of 
the standard cold dark matter scenario in an Einstein-de Sitter 
Universe from White et al. 1993), (ii) we limit the parameters 
that describe the geometry of the universe assuming that the 
gas fraction is constant in time, as firstly suggested by Sasaki 
(1996). 

The outline of our work is the following. In Section 2, we 
describe the cosmological framework that allows us to formu- 
late the cosmological dependence of the cluster gas mass frac- 
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tion. In Section 3, we use a sample of nearby clusters to con- 
strain mainly the cosmic matter density. Eight galaxy clusters 
with z > 0.7 are then presented and analyzed in Section 4 and 
a further constraint on the geometry of the Universe is given 
under the assumption of a constant gas fraction as function of 
redshift. A description of the systematic uncertainties that af- 
fect our estimates is discussed in Section 5. Finally, the com- 
bined probability function and the overall cosmological con- 
straints (also considered in combination with results from SN 
type la magnitude-redshift diagram) are described in Section 6 
along with prospective for future work. 

2. THE COSMOLOGICAL FRAMEWORK 

We refer to £l m as the total matter density (i.e., the sum of 
the cold and baryonic component: ft m — il c + f^) in unity 
of the critical density, p c = 3Hq /(8ttG), where H = lOOh 
km s _1 Mpc -1 is the Hubble constant and G is the gravita- 
tional constant, and to J7a as the constant energy density as- 
sociated with the "vacuum" (Carroll et al. 1992). We consider 
a generalization from this static, homogeneous energy compo- 
nent to a dynamical, spatially inhomogeneous form of energy 
with still a negative pressure, or "quintessence" (e.g. Turner & 
White 1997, Caldwell et al. 1998). We neglect the energy asso- 
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Fig. 1. (Left) Values of the density contrast A (with respect to what expected for an Einstein -de Sitter model) of the density 
contrast A estimated in an "f2 m + S7a = 1" universe at redshift 0, 0.1, 0.4, 0.7 and 1 (filled circles, from bottom to top, 
respectively) and with w equal to -0.2, -0.6, -1 (diamonds, from top to bottom, respectively; solid line: values at z = for 
w = — 1). (Right) Plot of the multiplicative factors that correct the gas fraction with respect to the initial estimate computed 

3 /2 

within a density contrast of 1500 in an Einstein -de Sitter universe: corrections due to the dependence upon d a n g (filled circles, at 
z =0.1, 0.4, 0.7 and 1, starting from the bottom, respectively) and the change due to the density contrast (diamonds; solid line: 
values at z = 0.1). 



ciated to the radiation of the cosmic microwave background, 
fi r w 5 x 10 -5 , and any possible contributions from light 
neutrinos, £1^ = n^/h 2 /92.5eV), that is expected to be 
less than 0.05 for a total mass in neutrinos, \i v , lower than 
2.5 eV (see, e.g., the recent constraints from 2dF galaxy sur- 
vey in Elgaroy et al. 2002 and from combined analysis of cos- 
mological datasets in Hannestad 2002). Thus, we can write 
O m + + f2k = lj where £lk accounts for the curvature 
of space. 

In this cosmological scenario, the angular diameter distance 
can be written as (e.g. Carroll, Press & Turner 1992, cf. eqn. 25) 

A — c S(uj) 
" an S ~~ H (l+z) |f2 k |i/2 i 

where S(u>) is sinh(cj), uj, sin(cj) for i\ greater than, equal to 
and less than 0, respectively, and 

e(z) = [n m (i + z f + n k (i + zf + n A (i + z) 3+3w ] 1/2 p.) 

that includes the dependence upon the ratio w between the 
pressure and the energy density in the equation of state of the 
dark energy component (Caldwell, Dave & Steinhardt 1998, 
Wang & Steinhardt 1998). Hereafter we consider a pressure- 
to-density ratio w constant in time (see, e.g., Huterer & Turner 
2001, Gerke & Efstathiou 2002 for the extension of eqn. [2] to 
a redshift-dependent form). In particular, the case for a cosmo- 
logical constant A requires w = — 1. 



2. 1. The cosmological dependence of the observed 
Cluster Gas Mass Fraction 

We assume that galaxy clusters are spherically symmetric grav- 
itationally bound systems. For each galaxy cluster observed at 
redshift z, we evaluate the gas mass fraction at r&, /gas (?*a) = 
Mgas{< r&)/M to t(< ta), where r& is defined according to 
the dark matter profile, M to t (< r), for a fixed density contrast 
A = M to t(< t A)/(47rp CjZ r'^). In the latter equation, p c z is 
the critical density at redshift z and is equal to / (8irG) 
with H z = HqE(z) (see eqn. |2j. In the following sections, 
we describe how the dark matter mass profile is obtained for 
each object and which density contrast we adopt initially in an 
Einstein-de Sitter universe. 

The assumed cosmological model affects the definition of 
the gas mass fraction, / ga s(^A)> given above in two indepen- 
dent ways: 



1. for a galaxy cluster observed at redshift z up to a char- 
acteristic angular radius 9 C and with an X-ray flux Sx = 
L x (l + z)-y{4ndl ng ) oc M^ s 9- 3 d- 3 g /dl ng [where 
the X-ray luminosity L x ~ n gas A(T gas ) x ^^ang> an( ^ 
A(T gas ) is the cooling function of the X-ray emitting 
plasma that depends only on the plasma temperature] and 
a total mass, M to t, estimated through the equation of the 
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hydrostatic equilibrium, the measured gas mass fraction is 

, _ M gas t> x V c flang , , q Ha,,/) 3 / 2 
■""tot Pc«ang 

(3) 

2. the density contrast, A, depends upon the redshift and the 
cosmological parameters. 

We have initially evaluated the gas fraction, / gas (As] m= i), 
in an Einstein-de Sitter universe with a Hubble constant of 50 
km s _1 Mpc -1 . Then, we change the set of cosmological pa- 
rameters and evaluate for each cluster at redshift z the new 
values of the angular diameter distance e? ang and the density 
contrast with respect to the critical density at that redshift. This 
density contrast is calculated according to the formula in Lokas 
& Hoffman (2001, Sect. 4.2) for a "fl m + fi A + O fe = 1" uni- 
verse and in Wang & Steinhardt (1998, eqn. 5, 7 and All; note 
that these formula are estimated for a background density that 
is Sl m (z) times the critical density) for a "quintessence" fiat 
model (see left panel in Fig. [3- 

Since we want to consider the gas fraction in each 
cluster estimated within the same A for any given set of 
cosmological parameters, we therefore multiply / ga s(Asi m =i) 
by two factors, the first one that rescales the distance, 

Fi = (dang,o m ^ A ,w/dang,f2 m =i) 3/2 ' and the second one 
that corrects by the change in the density contrast, F2 = 

(An m ,o A ,tu/Ao m =i) x (-fffd 2 ulg )n m =i/(-Hfd an g)n m ,o A ,w» 
where in the latter factor M to t(< ?"a) and the angular radius 
Qa = ^A/^ang are given (see right panel in Fig. Q. In 
particular, being A oc r^ 2 and / gas oc r^ 2 (e.g. Ettori & 
Fabian 1999a, and figure 13 in Frenk et al. 1999), we conclude 
that 

/gas,n m ,n A ,w — /gas,n m =i x Fi x F 2 . (4) 

As shown in Fig.Q the second correction affects the / gas 
values by less than 10 per cent and is marginal with respect to 
the cosmological effects due to the dependence upon the angu- 
lar diameter distance. We apply both these corrections to eval- 
uate each cluster gas fraction in the following analysis. 

3. FIRST CONSTRAINT: ft m FROM THE GAS 
FRACTION VALUE 

In this section, we describe how the local estimate of gas mass 
fraction provides a robust constraint on the cosmic matter den- 
sity. We make use of the further assumption of a constant gas 
fraction with redshift in the next section, where we consider a 
sample of galaxy clusters with z > 0.7. 

The observational constraints on the abundance of the light 
elements (e.g. D, 3 He, 4 He, 7 Li) in the scenario of the primor- 
dial nucleosynthesis gives a direct measurement of the baryon 
density with respect to the critical value, Sib- Moreover, the 
BOOMERANG, MAXIMA- 1 and DASI experiments have re- 
cently shown that the second peak in the angular power spec- 
trum of the cosmic microwave background provides a con- 
straint on Sib completely consistent with the one obtained 
from calculations on the primordial nucleosynthesis (e.g. de 
Bernardis et al. 2002). 



Table 1. The local sample from BeppoSAX MECS observa- 
tions. The quoted values are obtained from the deprojection 
of the spectral results and assuming a functional form of the 
total mass profile (see Ettori et al. 2002 for details). A Hubble 
constant of 50 km s _1 Mpc -1 is considered in an Einstein-de 
Sitter universe. 
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If the regions that collapse to form rich clusters maintain 
the same ratio f2b/O c as the rest of the Universe, a measure- 
ment of the cluster baryon fraction and an estimate of Sib can 
then be used to constraint the "cold", and more relaxed, compo- 
nent of the total matter density. This method alone can not pro- 
vide a reliable limit on the amount of the mass-energy presents 
in the Universe as hot constituents (e.g. WIMPS, like massive 
neutrino) or energy of the field (e.g. S1a, quintessence), as both 
do not cluster on scales below 50 Mpc. 

X-ray observations show that the dominant component of 
the luminous baryons is the X-ray emitting gas that falls into 
the cluster dark matter halo. Therefore, the gas fraction alone 
provides a reasonable upper limit on Sl c : 

Sic < oc h-^'\ (5) 

Jgas 

where the dependence of the ratio Slb// gas on the Hubble 
constant is factored out (White et al. 1993, White & Fabian 
1995, David, Jones & Forman 1995, Evrard 1997, Ettori & 
Fabian 1999a, Mohr, Mathiesen & Evrard 1999, Roussel, Sadat 
& Blanchard 2000, Erdogdu, Ettori & Lahav 2002, Allen, 
Schmidt & Fabian 2002). 

To assess this limit, we use the gas mass fraction esti- 
mated in nearby massive galaxy clusters selected to be relaxed, 
cooling-flow systems with mass-weighted T gas >4 keV from 
the sample presented in Ettori, De Grandi & Molendi (2002; 
cf. Table To date, this sample is the largest for which the 
physical quantities (i.e. profiles of gas density, temperature, 
luminosity, total mass, etc.) have all been derived simultane- 
ously from spatially-resolved spectroscopy of the same dataset 



4 



S. Ettori et al.: Cosmological constraints from the Cluster Gas Fraction 




0.05 0.10 0.15 0.20 
f gos (H =50 km s"' Mpc" 1 ) 

Fig. 2. Probability distributions (see Press 1996) at different 
density contrast A of the values of the gas mass fraction 
for cooling-flow clusters with T gas > 4 keV from the sam- 
ple in Ettori, De Grandi & Molendi (2002). The central and 
1 - a values are (0.126,0.024) and (0.152,0.020) at A = 
1500 and 500, respectively. These values have to be compared 
with the mean and standard deviation of (0.134,0.025) and 
(0. 152, 0.043), respectively. The estimated il^ from primordial 
nucleosynthesis results (Buries et al. 2001) is overplotted for 
comparison. 



(BeppoSAX observations, in this case). Through the deprojec- 
tion of the spectral results, and assuming a functional form for 
the dark matter distribution to be either a King (King 1962) 
or a Navarro, Frenk & White (1997) profile, the gas and total 
mass profiles are recovered in a self-consistent way. Hence, the 
density contrast, A = M to t(< T'A)/(47rp c . z r^), and the gas 
fraction at r A , /gasO^) = M gas (< r A )/M to t(< r A ), can be 
properly evaluated. 

In Fig. |2] we compare the estimated i7b from primordial nu- 
cleosynthesis calculation with the probability distribution (ob- 
tained following a Bayesian approach discussed in Press 1996) 
of the values of the gas mass fraction for nine (at A = 1500) 
and eight (at A = 500) nearby massive objects. This plot shows 
that at lower density contrast (i.e. larger radius) the amount 
of gas mass tends to increase relatively to the underlying dark 
matter distribution. Even if there is indication that the gas frac- 
tion becomes larger moving outward and within the virialized 
region, we decide to adopt / gas at A =500 as representative of 
the cluster gas fraction. In fact, at A =500 the gas fraction is 
expected to be not more than 10 per cent less than the universal 
value, a difference which is smaller than our statistical error, 
and in any case is likely to be swamped by other effects (see 
comments in the first item in Section 5). 



A = 500 



0.30 




o.oo r ..I i i 1 

40 50 60 70 80 90 
H 

Fig. 3. The observed gas fraction at A = 500 for the eight 
galaxy clusters in Table 1 is here plotted as function of the 
Hubble constant (shaded region: constraints from Freedman 
et al. 2001). The solid line indicates the central value and the 
dashed lines the 1 and 2 a uncertainties. 

We plot in Fig.[3]the dependence upon the Hubble constant 
of the / gas value (oc h~ 3 / 2 ) estimated at A = 500 for the ob- 
jects in the local sample. Assuming h — /ihst = 0.72 ± 0.08 
(from the results of the HST Key Project on distances mea- 
sured using Cepheid variables, Freedman et al. 2001) and 
il^h 2 = Slb.pN^ 2 = 0.019 ± 0.001 (from primeval deuterium 
abundance and calculations on the primordial nucleosynthesis, 
Buries, Nollett, Turner 2001), we obtain that fi c in eqn. |5]is 
less than 0.54 (95.4 per cent confidence level). 

Including a contribution from stars in galaxies of about 
/gal = 0.02(±0.01)/i 50 1 (White et al. 1993, Fukugita et al. 
1 998) and excluding any further components to the baryon bud- 
get (see e.g. Ettori 2001), one can write 51b/O c = /gas+/gal = 
/b and, consequently from the definition of fl m , f2b/^ m = 
/b/(l + /b). (Note that the estimate of M to t does not include 
the contribution of the gas mass, that would require the solution 
of a second order differential equation instead of a much sim- 
pler, and usually adopted, first order equation of the spherical 
hydrostatic equilibrium). 

Finally, we consider the eight relaxed nearby clusters i with 
T > 4 keV (see Tabled to evaluate the baryon fraction at red- 
shift Zi, fb,i, and within A(O m = 1, A = 0, w = -1) =500. 
For a given set of parameters (ft m , w) in the range [0, 1], 
[0,2] and [—1,0], respectively, we estimate /b,i (and its rela- 
tive error £b,i as propagation of the estimated error on / gas ,i 
and /gai, where the error on / gas ,i comes from the measured 
uncertainties on the gas and total mass estimates in Ettori et 
al. 2002, and the error on / ga i is Q-Olh^) after considering 
the cosmological dependence of both d ang and A(f2 m , w) 
(more relevant for high— z systems, see Section 2) and calcu- 
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Fig. 4. (Left) Probability distribution contours (solid lines: 1, 2, 3er for two interesting parameters) in the fib — Q m plane from 
marginalization of the likelihood provided from the baryon fraction in clusters assuming a /gal = 0.02(±0.01)/i5 1 (White etal. 
1993; Fukugita, Hogan & Peebles 1998), fi b ,p N (Buries et al. 2001), and H = 72 ± 8 km s" 1 Mpc" 1 (Freedman et al. 2001). 
As reference, fib,CMB (horizontal shaded region) and fi mj cMB (vertical shaded region; Netterfield et al. 2002), and (fib / fi m )2dF 
(dashed line indicates the central value; Percival et al. 2001) are indicated. (Right) Maximum likelihood distribution in the fi m 
axis after marginalization over the other parameters (dotted line: la, dashed line: 2a). 
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= E 

(O b -fl 



( /b.i-Ob/Om ) 2 



(/i-7ihst) 2 
71 ; 



(6) 



where fib.PN and /ihst are defined above. 

The x 2 -distribution in eqn. [6] is used to construct a A% 2 
statistics, A\ 2 — x\ ~ Xa min' hy which we generate regions 
and intervals of confidence (e.g., the la level of confidence for 
one and two degrees of freedom is A\ 2 = 1 and 2.3, respec- 
tively. A Xa ,min = 2 9-3 is obtained). 

Marginalizing over the accepted ranges of Hubble con- 
stant from the HST Key Project (Freedman et al. 2001) and 
fib from primordial nucleosynthesis, we obtain (2a) fl m = 



o.37i£:o8 and fi b 



0.032 



-0.017 

-o.oio 



, that are well in agreement 



with O m! cMB from CMB and fib/fi m from large scale struc- 
tures analysis of the "2dF" data (see Figure^. 

4. SECOND CONSTRAINT: GAS FRACTION 
CONSTANT IN TIME 

In this section, we present a sample of galaxy clusters with 
z > 0.7 and compute their gas mass fraction. These values 
are compared to the local estimates to put cosmological con- 
straints under the assumption that the gas fraction remains con- 
stant with redshift when computed at the same density contrast. 
It is worth noticing that this method, originally proposed from 



Sasaki (1996; see also Cooray 1998, Danos & Ue-Li Pen 1998, 
Rines et al. 1999, Ettori & Fabian 1999b, Allen et al. 2002), 
does not require any prior on the values of fib and Ho, tak- 
ing into account just the relative variation of the gas fraction as 
function of time. In other words, the method assumes that gas 
fraction in galaxy clusters can be used like a "standard candle" 
to measure the geometry of the Universe. This is a reasonable 
assumption in any hierarchical clustering scenario when the en- 
ergy of the ICM is dominated by the gravitational heating and 
is supported by numerical and semianalytical models for the 
thermodynamics of the ICM, also when including preheating 
and cooling effects. In recent hydrodynamical simulations with 
an entropy level of 50 keV cm" 2 generated in the cluster at 
z = 3 (see Borgani et al. 2002), the baryon fraction in clusters 
with an observed temperature around 3 keV is constant in time 
within few percent (see also figure 13 in Bialek, Evrard and 
Mohr 2001). Similar results are obtained in the semianalytical 
model of Tozzi & Norman (2001), where a constant entropy 
floor is initially present in the cosmic baryons. In Figure|5j we 
show the prediction for the baryonic fraction (in terms of the 
universal value) within an average overdensity A = 1500 as a 
function of the emission weighted temperature T ew , computed 
in a ACDM cosmology and with an entropy level of 0.3 x 10 34 
erg cm 2 g~ 5 / 3 (see Tozzi & Norman 2001 for details). A sig- 
nificant decrease of the baryonic fraction is expected at z = 1 
(dashed line) with respect to the local value (solid line) only for 
temperatures below 4 keV. In particular, at temperatures above 
6 keV, the baryonic fraction is constant or possibly 5 per cent 
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Fig. 5. The gas fraction in unit of the cosmic baryon bud- 
get at A = 1500 as a function of the observed (emission 
weighted) temperature computed for the model of Tozzi & 
Norman (2001) with a constant entropy of 0.3 x 10 34 erg cm 2 
g~ 5 / 3 in a ACDM (n m = 1 - fl A = 0.3) Universe. The solid 
line is for z — 0, and the dashed line for z — 1. Note the 20 
per cent offset from the universal value, which reduces to 10 
per cent for A = 500 (see comments on the baryonic depletion 
in Section 5). 

higher at z = 1 with respect to z = 0. This strongly supports 
the assumption of a baryon fraction constant with the cosmic 
epoch for clusters with kT > 4 keV. 

With this assumption, we expect to measure a constant av- 
erage gas fraction locally and in distant clusters. However, the 
gas fraction is given by a combination of the observed flux and 
of the angular distance, and thus it depends on cosmology (see 
discussion in Section 2). As shown in Figure ^ the high red- 
shift objects are more affected from this dependence and show 
lower / gas with respect to the local values when universes with 
high matter density are assumed. By requiring the measured 
gas fractions to be constant as a function of redshift, one can 
constrain the range of values of cosmological parameters which 
satisfies such a condition. 

4. 1. The high-z sample 

We define a local sample considering all the relaxed systems 
in Ettori et al. (2002) that have a mass-weighted temperature 
larger than 4 keV within A = 1500 (cf. Table [TJ. Within this 
density contrast, we measure / gas as described in Section 2. 
The density contrast is chosen to be 1500 as a good com- 
promize between the cluster regions directly observed in the 
nearby systems and those with relevant X-ray emission in 
the high— z objects. This sample includes eight high redshift 
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Fig. 6. Distribution as function of redshift of the gas fraction 
relative to the mean local value and estimated for the clusters in 
our sample at the same overdensity A =1500. Filled circles are 
the values calculated in an Einstein-de Sitter universe, whereas 
diamonds indicate the results for a low density universe (f2 m = 
1 - fi A = 0.3). 



(z > 0.7) hot (T > 4 keV) clusters, four of which selected 
from the ROSAT Deep Cluster Survey (RDCS; Rosati et al. 
1998, Stanford et al. 2001, Holden et al. 2002), two from the 
Einstein Extended Medium Sensitivity Survey (MS; Gioia et 
al. 1990), one from the Wide Angle ROSAT Pointed Survey 
(WARPS; Perlman et al. 2002, Maughan et al. 2002) and one 
part of the North Ecliptic Pole survey (NEP; Gioia et al. 1999, 
Henry et al. 2001). We reprocess the level=l events files re- 
trieved from the archive and obtain a spectrum and an image for 
each cluster (see details in Tozzi et al., in preparation). Seven of 
these objects were observed in ACIS-I mode (only MS 1054 has 
been observed with the back-illuminated S3 CCD). Following 
the prescription in Markevitch & Vikhlinin (2001), the effec- 
tive area below 1.8 keV in the front-illuminated CCD is cor- 
rected by a factor 0.93 to improve the cross calibration with 
back-illuminated CCDs. The spectrum extracted up to the ra- 
dius r out to optimize the signal-to-noise ratio is modelled be- 
tween 0.8 and 7 keV with an absorbed optically-thin plasma 
(wabs (mekal) inXSPECv. 1 1.1.0, Arnaud 1996) with fixed 
redshift, galactic absorption (from radio HI maps in Dickey & 
Lockman 1990) and metallicity (0.3 times the solar values in 
Anders & Grevesse 1989) and using a local background ob- 
tained from regions of the same CCD free of any point source. 
The gas temperature and the normalization K of the thermal 
component are the only free parameters. The surface bright- 
ness profile obtained from the image is fitted with an isother- 
mal j3— model (Cavaliere & Fusco-Femiano 1976, Ettori 2000), 
which provides an analytic expression for the gas density and 
total mass of the cluster. In particular, the central electron den- 
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Table 2. The high redshift sample from Chandra observations. The results are obtained applying a 0— model in an Einstein-de 
Sitter universe with a Hubble constant of 50 km s _1 Mpc -1 . All the quoted errors are at la level. Note that MS1054. 5-0321 
presents significant substructure (e.g. Jeltema et al. 2001). The temperature and the best-fit of the surface brightness profile are 
estimated from the main body of the cluster once a circular region centered at (RA, Dec; 2000)=(10 h 56 m 55"7, — 3°37'37") and 
with radius of 36 arcsec is masked. 
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sity is obtained from the combination of the best-fit results from 
the spectral and imaging analyses as follows: 



2 47Td2 ng X (1 + zf X K X 10 14 

n 0,ele = ~~ 



0.82 x 4rrr3 x EI 



(7) 



where the emission integral is estimated by integrating along 
the line of sight the emission from the spherical source up to 
lOMpc, EI= f* 1 (l + x 2 )- 3 f 3 x 2 dx+ f* 2 (l + x 2 )- 3 ^ 2 (l- 
cos6)dx, with 9 = arcsin(r out /r), x\ = r out /r c and X2 — 
10Mpc/r c , (/?, r c ) are the best-fit parameters of the (3— model 
and we assume n p = 0.82n o in the ionized intra-cluster 
plasma. After 1000 random selections of a temperature, nor- 
malization and surface brightness profile (drawn from Gaussian 
distributions with mean and variance in accordance to the best- 
fit results), we obtain a distribution of the estimates of the gas 
and total mass and of the gas mass fraction. We adopt for each 
cluster the median value and the 16th and 84th percentile as 
central and la value, respectively (see Table|2j- 

In Figure [6J we plot / gas of the clusters in exam as a func- 
tion of redshift, as computed in an Einstein-de Sitter (filled 
dots) and a low density (Q m = 1 — Oa = 0.3) universe (open 
squares). 



4.2. The analysis 

We compare our local estimate of the gas mass fraction with 
the values observed in the objects at z > 0.7 within the same 
density contrast A = M tot (< r A )/ (47rp c ^ z r^), where M to t(< 
ta) is estimated from the (3— model. Initially, we estimate the 
gas fraction for all the clusters at A = 1500 in an Einstein-de 
Sitter universe. Then, we proceed as described in Section 2. 

As discussed above, we consider a set of parameters 
(fi m , Oa) in the range [0, 1] and [0, 2], respectively, both fix- 
ing w equals to —1 as prescribed for the "cosmological con- 
stant" case and exploring the range w G [—1,0] (in this case, 
we require flk = 0), and evaluate the distribution 



2 \ ^ I /gasj (^mi S^a) /g 

xb = }^ j — , ^-2 



(8) 



gas,j 



c gas 



where / gas and are the mean and the standard deviation of 
the values of the gas fraction in the local cluster sample, and 
£gas,j is the error on the measurement of / gas ,j for j £ [high-z 
sample] . It is worth noticing that the use of the standard devi- 
ation around the mean is a conservative approach. For exam- 
ple, at (A, fi m , Hq) = (1500,1,50), we measure a mean of 
0.134 and standard deviation of 0.025, the latter being about 8 
times larger than the measured error on the weighted mean and 
slightly larger than the dispersion obtained with the Bayesian 
method illustrated in Fig. 12 Using just the results obtained on 
the Xb distribution, we obtain a best-fit solution that requires 
the following upper limits at 2a level (one interesting parame- 
ter), fi m < 0.64 and Q A < 1.69 (cf. Figure|7}. 



5. SYSTEMATIC UNCERTAINTIES 

The intrinsic scatter in the distribution of the / gas values of 
approximately 20 per cent is the most relevant statistical uncer- 
tainty affecting our cosmological constraints. We discuss in this 
section how a number of systematic uncertainties contribute to 
a lesser extent. 

- We have assumed that the intracluster medium is in hy- 
drostatic equilibrium, distributed with a spherical geome- 
try, with no significant (i) dumpiness in the X-ray emitting 
plasma, (ii) depletion of the cosmic baryon budget at the 
reference radius, and (iii) contribution from non-thermal 
components. The hydrostatic equilibrium in a spherical po- 
tential is widely verified to be a correct assumption for local 
clusters, but it cannot be the case (in particular on the geom- 
etry of the plasma distribution, but see Buote & Canizares 
1996, Piffaretti, Jetzer & Schindler 2002) for high red- 
shift clusters. The level of dumpiness in the plasma ex- 
pected from numerical simulations (e.g. Mohr et al. 1999) 
induce an overestimate of the gas fraction. On the other 
hand, the expected baryonic depletion (still from simula- 
tions; e.g. Frenk et al. 1999) underestimates the observed 
cluster baryon budget by an amount that can be compara- 
ble (but of opposite sign) to the effect of the dumpiness 
(see discussion in Ettori 2001). For the high redshift ob- 
jects observed within A=1500, a larger dumpiness can be 
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Fig. 7. Maximum likelihood distributions in the "fi m + + 
Qf. = 1" region obtained from the application of the method 
that requires "/ gas =constant" (cf. Section 4; shaded region: no- 
big-bang solution, dotted line: "0^ = 0" region). The con- 
tours enclose the regions with A\ 2 =2.30, 6.17, 11.8 (with 
respect to the minimum of 11.7 with eight objects in exam), 
corresponding to 1, 2, 3 er, respectively, for a distribution with 
two interesting parameters. 



present as consequence of on-going process of formation 
that, on the other hand, might still partially compensate for 
a residual depletions of baryons observed in the simulations 
in the range 4-6 keV (cfr. Figure |5Jl. The combination of 
these effects induce however an uncertainty in the relative 
amount of baryons of less than 10 per cent, percentage that 
is below the observed statistical uncertainties. 
Moreover, the radial dependence of / gas can be slightly 
steeper in the inner cluster regions due to the relative 
broader distribution of the gas with respect to the dark mat- 
ter (but see Allen et al. 2002 for / gas profiles observed 
flat within A=2500). Finally, from the remarkable match 
between the gravitational mass profiles obtained indepen- 
dently from X-ray and lensing analysis (e.g. Allen, Schmidt 
& Fabian 2001), it seems negligible any non-thermal con- 
tribution to the total mass apart from a possible role played 
in the inner regions. 
- The high redshift sample has been analyzed using an 
isothermal (3— model. From recent analyses of nearby clus- 
ters with spatially resolved temperature profiles, it appears 
that this method can still provide a reasonable constraint on 
the gas density, but surely affect the estimate of the total 
mass, in particular in the outskirts and, more dramatically, 
when an extrapolation is performed. On the other hand, 
the physical properties of nearby clusters at overdensity of 
1500 (e.g. Ettori et al. 2002) guarantee that the expected 
gradient in temperature is not significantly steep and, there- 



fore, does not affect the estimate of the total gravitating 
mass. However, the presence of a gradient in the tempera- 
ture profile would reduce the total mass measurements and 
increase, consequently, the derived gas mass fraction. 
- We assume that relaxed nearby and high-z clusters, both 
with T gas > 4 keV, represent a homogeneous class of ob- 
jects. If we relax this assumption and include non-relaxed 
nearby systems (or equivalently no-cooling-flow clusters in 
the sample of Ettori et al. 2002), we obtain the same cos- 
mological constraints though the scatter in the distribution 
of the gas fraction values becomes larger (e.g., at A =1500, 
the scatter is 0.025 for relaxed systems only and 0.033 for 
the complete sample). 

In conclusion, we check our results against two effects: (i) 
by increasing the gas mass fraction in the high— z sample by 
a factor 1.15, and (ii) by changing the slope of the radial de- 
pendence of /gastj) (cfr. factor Fi in Section 2) by ±0.1. All 
these corrections do not change significantly the results quoted 
in the next section on O m due to heavier statistical weight given 
from the first method, that is the one less (or completely not) 
affected from the above mentioned effects. On the other hand, 
they mostly affect the second method in such a way that (i) gas 
fraction values at high— z higher by 15 per cent reduce the 2a 
upper limit on 57a by 10 per cent, but increase the one on w 
by 50 per cent, (ii) steeper radial / gas profiles (i.e. larger de- 
pendence upon the density contrast A) raise the limit on 
but have not relevant effects on w. In details, when the slope 
changes from to —0.1 and —0.2, the upper limit on in- 
creases from 25 to 5 per cent, respectively. 

6. CONCLUSIONS 

We show how the combined likelihood analysis of (i) the rep- 
resentative value of /g as in clusters of galaxies and (ii) the 
requirement that f gas {z) =constant for an assumed homoge- 
neous class of objects with T > 4 keV can set stringent limits 
on the dark matter density and any further contribution to the 
cosmic energy, i.e. J7 m and tt\ respectively. 

First, a total x 2 distribution is obtained by combining the 
two distributions presented in eqn. |6]and |8] i.e. \ 2 = x\ + 
Xb- The resulting likelihood contours (Figure [8} are obtained 
marginalizing over the range of parameters not investigated. 
With further a priori assumptions on fib and Ho, and assuming 
a flat geometry of the universe, we constrain (see right panel in 
Figure[8j the dark energy pressure-to-density ratio to be 

w < -0.82(lcr), -0.49(2cr), -0.17(3ct). (9) 

This constraint is in excellent agreement with the bound on w 
obtained with independent cosmological datasets, such as the 
angular power spectra of the Cosmic Microwave Background 
(e.g. Baccigalupi et al. 2002), the magnitude-redshift relation 
probed by distant type la Supernovae and the power spectrum 
obtained from the galaxy distribution in the two-degrees-field 
(for a combined analysis of these datasets, see, e.g., Hannestad 
& Mortsell 2002 and reference therein). Moreover, this upper 
bound is completely in agreement with w = —1 as required 
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Fig. 8. (Left) Maximum likelihood distributions in the "fi m + + &k = 1" region. Contour plots (thick solid lines) from the 
combination of the two likelihood distributions (A: cluster baryonic content, B: gas fraction constant with redshift; dashed lines 
indicate the constraints from the second method only, see Fig.0 the cross indicates the best-fit result at f2 m = 0.34, £1 A = 1.30) 
with overplotted the constraints from the magnitude-redshift method applied to a set of SN la (cf. Leibundgut 2001, thin solid 
lines). (Right) Constraints on the parameter w of the cosmological equation of state (thick solid lines from the combination of the 
method A and B; dashed lines from the second method only; the crosses show the best-fit results, that is located at Q m = 0.35 
and w — — 1 for the combined probability distribution). The thin solid lines indicate the constraints from SN la (P. Garnavich, 
priv. comm.; updated version of Garnavich et al. 1998 combining Riess et al. 1998 and Perlmutter et al. 1999 dat sets). The 
contours enclose the regions with A% 2 =2.30, 6.17, 11.8, corresponding to 1, 2, 3 a, respectively, for a distribution with two 
degrees-of-freedom. 



for the equation of state of the "cosmological constant". Fixing 
w = — 1, we obtain (for one interesting parameter) 

Sl m = 0.34t°;°3(l«7), ±° f 5 (2a), ±°;&(3<r) 

A = 1.30t°j°(la), ±?;^(2a) ) +^(3a). (10) 

Finally, imposing a flat Universe (i.e. = 0) as the recent 
constraints from the angular power spectrum of the cosmic mi- 
crowave background indicate (e.g. de Bernardis et al. 2002 and 
references therein), we obtain that Q m = 1 — Q\ = 0.33^oo5 
at 95.4 per cent confidence level. 

Our limits on cosmological parameters fit nicely in the cos- 
mic concordance scenario (Bahcall et al. 1999, Wang et al. 
2000), with a remarkable good agreement with independent es- 
timates derived from the angular power spectrum of Cosmic 
Microwave Background (Netterfield et al. 2002, Sievers et al. 
2002), the magnitude-redshift relation for distant supernovae 
type la (Riess et al. 1998, Perlmutter et al. 1999; the like- 
lihood region from a sample of SN type la as described in 
Leibundgut 2001 is shown in Figure [8] panel at the bottom- 
left), the power spectrum from the galaxy distribution in the 
2dF Galaxy Redshift Survey (e.g., Efstathiou et al. 2002) and 
from galaxy clusters (e.g. Schuecker et al. 2002), the evolution 
of the X-ray properties of clusters of galaxies (e.g. Borgani et 
al. 2001, Arnaud, Aghanim & Neumann 2002, Henry 2002, 
Rosati, Borgani & Norman 2002). For example, combining 



the constraints in Figure [8] between the allowed regions from 
the gas mass fraction and the magnitude-redshift relation for 
SN-Ia, we obtain (2 a statistical error) il m = 0.34+QQg and 
51a = 0.94+Q32- These values are 0.5 and 0.6 er higher, re- 
spectively, than the CMB constraints obtained with the SN- 
Ia prior (see Table 4 in Netterfield et al. 2002). Moreover, by 
combining / gas and SN-Ia measurements we can obtain a very 
tight constraint on w (right panel in Figure [8): w < —0.89 and 
il m — 0.32+Q Q5 at the 95.4 per cent confidence level. 



We have demonstrated how the measurements of the cluster 
gas mass fraction represent a powerful tool to constrain the cos- 
mological parameters and, in particular, the cosmic matter den- 
sity, n m . Nonetheless, the limits on 17 a and w, though weaker, 
provide a complementary and independent estimate with re- 
spect to the most recent experiments in this field. On this item, 
it is worth noticing that our constraints on 57a are mostly due 
to the d ang dependence of / gas (cf. Fig.[0- Thus, a larger sam- 
ple of high— z clusters with accurate measurements of the gas 
mass fraction will significantly shrink the confidence contours, 
as we show in Figure [9] Compilations of such datasets will be 
possible in the near future using moderate-to-large area surveys 
obtained from observations with Chandra and XMM-Newton 
satellites. 



10 



S. Ettori et al.: Cosmological constraints from the Cluster Gas Fraction 



2.0 



< 
a 



0.5 



5 obj, z=[0.8,1.2] 20 obj, z=[1,2] 
5 obj, {=[1.2,1.6] 40 obj, z=[l,2] 




0.2 0.4 0.6 0.8 

n 



1.0 



Fig. 9. Constraints in the fi m — SI a plane as in Fig. for 
simulated sets of gas mass fraction measurements. The con- 
tours are at 1 a confidence level. A cosmological model with 

3 /2 

fl m = 1 — S2a = 0.34 is assumed and only the c? a ng depen- 
dence is considered. In the redshift range [1, 2], the increase by 
a factor of two of the sample allows to reduce the upper limit 
on J7a by about 10 per cent. 
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Fig. A-l. Data and best-fit MEKAL model of the spectrum of the galaxy clusters at high redshift in our sample. 




Fig. A-2. Data and best-fit /3-model of the surface brightness profile of the galaxy clusters at high redshift in our sample. Dotted 
lines indicate the best-fit background value. 



